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_ ABSTRACT 


An iterative technique for the retrieval of temperatures 
at each of 100 levels ranging from 1000 mb to 0.01 mb is 
evaluated. Clear-column radiance data in the carbon dioxide 
channels of the VTPR of NOAA 2 are used in inverting the 
radiative transfer equations to deduce the T(P) profile. MThe 
retrieval technique includes the computation of atmospheric 
transmittances due to three atmospheric absorber masses 
(carbon dioxide, water vapor, and ozone) and non-homogeneous 
temperature-pressure effects along the vertical. The program 
also corrects these transmittances for zenith path differences 
between the satellite and the retrieval site when the site is 
not directly below the Sensor. A normaiized wave number, 


- but decreasing with decreasing values of 


close to 700 cm 
pressure, was used in solving the Planck function for tempera- 
ture, level-by-level. The retrieved temperatures were then 
compared with matching radiosondes at 15 key atmospheric 


levels and the resulting errors analyzed on a globally averaged 


emadelatituae band basis at each of the 15 levels. 
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I. INTRODUCTION 


The quantitative measurement of temperatures at known 
pressure levels in the atmosphere is fundamental to present 
Gay weather analysis and forecasting. Until the launch of 
the Nimbus 3 satellite with its Satellite Infrared Spectrometer 
(SIRS) in April 1969, free air temperature soundings were 
obtained entirely from balloons, aircraft, and rockets. Since 
these sources were readily available only from populous land 
areas and a few ships, large areas of the world were without 
any reports. World-wide distribution of temperature sounding 
fee was mace possible by the use of satellites. The retrieval 
of a temperature profile by remote-sensing satellite measure- 
ments was first suggested by Kaplan (1959), and has been the 
object of numerous research studies. 

The Nimbus 3 and Nimbus 4 satellites were equipped with 
SIRS-A and SIRS-B respectively. These first generation SIRS 
instruments, SIRS-A with a vertical scan, and SIRS-B with an 
additional side scan, measured monochromatic radiances in seven 
narrow regions or channels of the 15 yum band of carbon dioxide, 
and one in the window region of water vapor centered at 
B99 cm +, Figure 1 (from Gelman et al. 1972) shows the 
weighting functions dt;/d(inP), after Smith et al. (1972), 
of the radiometric emitting channels of the SIRS-B instruments. 
A peak in dt, /dx Pomatizemto thie InP axis depicts the layer 


in the vertical contributing the maximum radiance to the 
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feeeellite channel sensor. These radiances could tnen. pe 


Beveted tO temperature profiles through the Radiative Transfer 
Equation (RTE) 


P 


Oo ole 
Be a 
I. = BIT (P.)] eee ey) + mi Bree) ) ap. aP : (1) 
= | 
Here, I. 1S the spectral radiance in channel 1 (i=1,2,...,8), 


1 
B,[T(P)] is the well-known Planck function for wave number Vs 


and at temperature T(P), T; (P) is the atmospheric transmittance 
averaged over wave number i (where i signifies the channel 


number centered at Vi), and P is the pressure-level for which 


tT; (P) has been computed. 
001 


| 668.7¢m' 


25 
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BIGURE 1. Derivative Of transmittance (1) with respect 
yy . 
to x(P) = P la for SIRS-B 


(after Gelman et al. 1972) 





Retrieval of temperature profiles from the SIRS-B radi- 
ometer was continued in systematic fashion by the National 
Environmental Satellite Service (NESS) until the radiometer 
output was considered non-usSable on 8 April 1971. The number 
of retrievable satellite soundings collected during 1969-71 
was insufficient to make a marked impact on the analysis 
problem over the oceans. 

Anticipating an increaSed demand for quantitative weather 
observations from orbiting satellites, the National Oceanic 
and Atmospheric Administration (NOAA) in conjunction with the 
National Aeronautics and Space Administration (NASA), developed 
a second generation set of remote-Sensing radiometers, improved 
over the SIRS-type, that would provide twice daily coverage 
of most of the globe and a higher resolution at individual 
mean Spots. 

As a result of these developments, the first of the new 
instruments was the Infrared Temperature Profile Radiometer 
(ITPR) as used on Nimbus 5. This radiometer carried four 
carbon dioxide channels, two window channels at 4.3 um and 
11.5 um, respectively, and a water vapor channel near 20 um. 

In 1972, NASA launched the NOAA-2 satellite, which carried 
aloft a new remote-sounding radiometer for which the ITPR 

had been a prototype: The new instrument, the Vertical 
Temperature Profile Radiometer (VTPR), had been Still further 
improved over the ITPR in its angular-array scanning technique. 


The VTPR carried six carbon dioxide channels, a window channel, 
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mca water vapor channel. Observations made at scan spot 
arrays in channels 1,2,...,8 made it possible to systematically 
ee ieve clear-column radiances at a higher resolution than 
before. The new instrument was designed to provide two tem- 
perature profiles per grid point per day on a 400 kilometer 
spacing over the oceans, or a global total of 6000 temperature 
profiles per day. An actual yield of 72%, or about 4000 
profiles per day has been realized [Jastrow and Halem (1973)]. 
The use of these VTPR clear-column radiances from the NOAA 2 
satellite and the temperature retrieval system described by 
McMillan et al. (1973) was the basis for a thesis by Moran 
(1974). The basic differences between the Moran study and 
previous work by Smith et al. (1972) was the use, by Moran, 
of standard atmospheric carbon dioxide transmittances with no 
interactive effects on the channel 1,2,...,6 transmittances 
from the other atmospheric constituents. Moran aiso computed 
radiances by the use of a 17-layer quadrature scheme for inte- 
Seatton Of Eq. 1 rather than by means of the more usual 100- 
level atmosphere. This was done in order to attempt to reduce 
the trunction error involved with uncertainties in the integra- 
fron Of the function B, ([T(P)]dt,/dp when a relatively crude 
model for Tt; (P) was utilized. 

The first objective of the present thesis was to improve 
upon the Moran retrieval technique by incorporating recently 
updated state-of-the-art methods [Fritz et al (1972)] for 


computing atmospheric transmittances, including the combined 
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effects of water vapor, ozone, and carbon dioxide in each 
channel, aS well as their vertical path dependencies on 

Pe ecure and temperature. It was also necesSary to correct 
these transmittances in each channel for the applicable zenith 
angles of the clear-column spot. For this purpose, a Special 
data set for March 1973, including the tuned transmittance 
coefficients for the period, coupled with the first-guess 
temperature and humidity profiles, was acquired from NESS. 

A second and equally important objective was to test the 
final temperature retrievals by comparison with radiosondes 
at island and coastal stations matching the VTPR scan spots 
in both time and spatial location. The retrieved profiles 
were compared at each of the 15 standard pressure levels 
(Table 3) with the matching radiosondes and were analyzed 
both by latitude and by pressure level. 

The retrieval system used here has the feature of inte- 
grating Eq. 1 over 33 atmospheric layers at each clear-column 
radiance scan spot. This system was incorporated in order 
to minimize the truncation error in the quadrature scheme 
applied to the function B, dt, /dp which appears in the RTE 
(Eq. 1). The improved accuracy in the determination of the 
temperature profile using this method is one of the main 


merusts of this thesis. 
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fi Saul TE DATA 


The NOAA series of satellites are modifications of the 
earlier Environmental Science Service Administration (ESSA) 
Satellites. The major change was the elimination of all 
vidicon camera systems in favor of radiation-sensing instru- 
ments, and in particular the inclusion of the vertical 


temperature profile radiometer (VTPR). 


Pee NOAA 2 DESCRIPTION 

The NOAA 2 satellite orbits the earth at an altitude of 
1464 kilometers once during each 115 minutes. Figure 2 
illustrates the projection of seven orbits onto the earth. 

The solid lines indicate the north-to-south or "descending" 
portions of the orbits, which occur over the sunlit hemisphere; 
dashed lines are nighttime "ascending" portions of the same 
Orbits. Equatorial crossings occur at 0900 aeecenaina) and 
2100 (ascending) local solar time. 

The VTPR instrument scans perpendicularly relative to the 
Subsatellite track from left to right, when viewed downtrack, 
in 23 individual non-overlapping steps covering +30.3 degrees 
from local nadir. The instrument observes each spot in eight 
spectral intervals or channels through the use of the automatic 
optical filter selector equipment of the VTPR. Six of the 
channels are located in the 15 um band of carbon dioxide, one 


in the 19 um region of the rotational band of water vapor, and 
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one in the 12 um atmospheric window (see Table 1). After 
the entire set of 8 spectral measurements has been obtained, 
the Scanning mirror rotates to the next spot. 

Table 1. VTPR channel designators and centroid wave numbers 


EOpemieleet set spour Of NOAA] 2 VTPR instrument 
momMoer Im iusead) in March 1973). 


CHANNEL i BD 3 4 5 6 7 8 

WAVE NO. 668.5 677.5 695.7 707.2 724.7 747.5 534.5 835.1 
-1 

(cm ~) 


B. CLEAR-COLUMN PROCEDURE 

In order to obtain an accurate temperature retrieval, 
the “raw" radiances, that is the radiances measured at the 
VTPR instrument, must be adjusted to eliminate the effects 
of any clouds that may be located in the field of view. 
Smith (1970) and others introduced cloud layer models wherein 
the extent and thickness of cloud layers which would give the 
meadiance aS Sensed at the Satellite are mathematically intro- 
Guced into the temperature profile calculation. However, the 
results in the neighborhood of a single scan spot tend to 
result in non-unigque cloud heights and amounts which causes 
uncertainty in the deduced clear-column radiances. The 
method described below was later developed by Smith (1971) 
in order to eliminate these uncertainties. 


ThesNational Environmental Satellite Service (NESS) 


eliminates cloud contamination through the use of a conversion 


program called CLRAD, applied to a large array of contiguous 
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raw radiance spots. Sets of eight scan lines, containing 
23 scan spots are divided into 3 subarrays or boxes of 8 by 
Epcos by 7, and 8 by 8 scan spots as shown in Figure 3. Ob- 
servations of raw VTPR radiances are made in each of the eight 
channels and at each scan spot. Figure 4 shows the lateral 
Scan geometry along each scan line, where spot 5, spot 12, 
and spot 19 in line 4, are the center of box left, box center 
and box right, respectively. Each box-center is denoted by 
an X in Figure 3. At these locations it is necesSary to 
determine the exact zenith angles to the sensor aboard the 
satellite using the following relationship 

SIN(Z) = SIN(N,) 3= (2) 
Here Z is the zenith angle at scan spot n (n = 5,12,19). 


N is the nadir angle at the same spots, a is the constant 


"6370 km (earth's radius), and H is the altitude of the 


satellite (mean value of 1460 km). The mean zenith angle at 
spot 12 is 0°, while at spots 5 and 19 it is 23°47' left and 
right of center respectively. 

In addition to the VTPR radiances, meaSurements in a 
Window-channel Scanning Radiometer (SR), sensing in 10.5 to 
12.5 um, permits determination of the sea surface temperatures 
at each of the scan spots. McMillin et al. (p 25-33, 1973) 
outline a process of extrapolation of raw-radiance to clear 


column radiance within localized sub-areas about the scan 
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Figure 3 — VTPR scan pattern and data analysis array 
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Figure 4, NOAA-2 VTPR lateral scan geometry at raw 
radiance spots 
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spots. The specifics of the clear radiance conversion method 
involves: (1) the specifications of a reasonably accurate 
value of the first guess temperature profile T(P) to 10 mb, 
(2) an estimation of the water vapor transmittances in the 
Window channel derived from the T(P), and (3) the Fee eee 
tion of the sea surface temperature derived from the SR. The 
value of the clear-column window channel radiance can then be 
computed and compared with the observed window channel radi- 
ance. If the observed “raw" radiance iS equal to or greater 
than the computed radiance, the sounding is considered to 
be cloud free. If the computed radiance is greater than the 
"raw" radiance, the sounding is considered to be cloud 
contaminated and the extrapolation process 1s necesSary. 

To facilitate the extraction of the clear-column radiance 
set, “raw" radiance values from adjacent spot pairs (spots 1 
and 2 with their associated T, (vy) and I, (v) for example) are 


extrapolated to I (v). This can be done if the cloud 


clear 
amounts vary linearly in the field of view and share approxi- 
mately the same cloud top level. Figure 5 shows the patterns 
used by NESS for pair comparisons. Each 8 by 8 box, so 
analyzed, will contain 49 scan spots, which can be compared 


with four adjacent spots for a total possible yield of 196 


(o) 2) Here 1 denotes 


values of clear-col aLsie crs) if 
clear-column radianc alco al 


all channels other than the window channel since the window 
channel gives rise only to a possible 49 clear-column esti- 
mates. A proportionately smaller number of both radiance 


sets are found in the 8 by 7 box. If a suitable number of 
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the scan spots are capable of giving. four adjacent clear- 
cOlumn values without excesSive variances, all resulting 
values are usable in obtaining the weighted clear-column 


radiance in the respective channels. 


STEP BY STEP 
ROPE LNs wi 


STEP BY STEP 
TO SPOT 22 


Figure 5. Scanning patterns used by the National Environmental 
Satellite Service (NESS) for pair comparison in 
determining VTPR clear-column radiance. 
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Research on refining the accuracy Srrcetermindtion Of 
clear-column radiances is actively being conducted. This 
involves both insStrumental accuracy improvements as well as 
the algorithms for the extrapolations of cloud amounts to 


Zero. 





ce es PrRAtiotlt TANGE MODELS IN THE VITPR CHANNELS 

A. ATMOSPHERIC TRANSMITTANCE COMPUTATIONS: GENERAL 

CONSIDERATIONS 

Atmospheric transmittances for the Six carbon dioxide 
Spectral intervals used in the retrieval are fundamental and 
must be computed accurately. The resultant transmittance in 
any given channel and at any pressure level P(J) in the atmos- 
phere may be expressed as the product of the transmittances 
Gue to the carbon dioxide alone Tt[C(J)] and the overlapping 
transmittances due to the integrated ozone, OZ(J), and water 
vapor, U(J), from the top of the atmosphere to the level J 


respectively. Thus 
Eo) Sa CCl (U0 (di) ]*7 [Oz (J) ] (3) 


Here C(J) denotes the total absorber mass of the carbon 
mu@arde from P(l) = 0.01 mb to P(J). The symbol Oz(J) ce- 
notes the corresponding absorber mass of ozone expressed in 
NTP Standard amounts per square centimeter column above level 
J, and U(J) is the precipitable water vapor in the same column. 
As an example of the transmittances produced by the program 
discussed in this thesis, sample outputs are included following 
the program 1isting. 

Moran (1974) used only carbon dioxide transmittances 
acting in the vertical column (i.e., {U(J) J=1 [02 (J) J=1, by 


assumption) and further simplified the computation by selecting 
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a mid-latitude, standard-atmosphere .transmittance model for 
carbon dioxide in the Six listed channels. Moran found that 
teer-adiabatic temperature lapSe rates tended to occur 
between the 1000 mb and 900 mb levels. ConSequently he found 
1t necesSary to adjust the values of the transmittance profile 
somewhat arbitrarily in channels 5 and 6, Since the trans- 
mittance from the surface in these channels tended to be 
somewhat higher than with the more realistic product 
transmittances of Eq. 3. 

Transmittances in the carbon dioxide channels were based 
upon the laboratory obServations of the carbon dioxide infra- 
red Spectrum and resultant calculations by Drayson (1966), 
made on a line-by-line basis asSuming a homogeneous path 
meor a fixed temperature). 

Mimeilewlattons applied tO nen-nomogeneeus vertical paths, 
Smith (1969), calculates the carbon dioxide transmittances 
corresponding to the U. S. Standard Atmosphere (1962); and 
makes improvements on the usual Curtis-Godson modifications 
[Goody (1964)] applied to path transmittances for varying 
presSure, temperature, and absorber masses. In addition, 
modified profiles corresponding to temperature differences 
meem the U.S. Standard Atmosphere by +410, +20, and +30°K 
were also computed by Smith, after Drayson's (1966) laboratory 
calculations. A method for applying temperature corrections 
to the abSorber maSs in each non-homogeneous sub-layer was 
devised so that the effective temperature-modified carbon 


dioxide transmittances could be adjusted by a variation of 
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the Curtis-Godson method applied to successive layers, uSing 
the laboratory data of Drayson. 

| The method utilized here for carbon dioxide temperature 
adjustment along the path will be discussed in greater detail 
following that for water vapor. This latter method is some- 
what more general due to the variable mixing ratio of this 
constituent in the vertical as well as the greater complexity 
of the individual line structure in the water vapor infrared 


Seeccrum. 


Be WATER VAPOR = MIXING RATIO ESTIMATION PROCEDURE 

In order to determine the transmissivity due to water 
vapor in each of the VTPR channels, the precipitable water 
vapor integrated from J = 01 to any arbitrary level must be 
determined. For this purpose an estimation procedure for 
determining the mixing ratio , W(J), at each level was used 
to represent the water vapor profile. The parameter W(J) 
was used because it was found to be determined from T(P) 
alone with reasonable skill, using the method of Weinreb 
femerosoy (1973); and also because it enters directly into 
the calculation of the water vapor transmittances in the 
eight VTPR channels. 

The calculation of the mixing ratio W(J) may be derived 


from 


W(J) = SAT(T,J) *H(J) (4) 
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H(J) is the relative humidity and is an unknown function of 


‘) 


pressure. The use of the water vapor channel (v = 535 cm. 
was not attempted for retrieval in this thesis. Furthermore, 
wort. og), the saturated mixing ratio is a function of 


temperature T(J) and presSure P(J) only, and is calculated 


by the GOFF-GRATCH formula [List 1963] 


log 6, = -7.90298(T,/T-1) + 5.02808 log(T_/T) (5) 
= 1.3816x10 /f1ott-34441-T/Ts) _ 4 
men lon 5 8’) =) 2) + jog ec 


WS 


Here e 1s the saturation vapor pressure over a plane water 
Surface, 


T is the temperature in °K, 
T is the steam point temperature 373.16 °K, and 


e 1s the saturation vapor presSure of water at the steam 
point temperature. 


The value of SAT(T,J) is then computed uSing the relationship 


e 
W 


Since no equation directly relates the relative humidity 
to the temperature profile, the mixing ratio was estimated 
from the saturated mixing ratios at a specified set of eleven 
levels (Table 2) by a least-squares formulation based solely 
upon these eleven saturated mixing ratio predictors after 


Weinreb and Crosby (1973). The method formulates the mixing 
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ratio W(J) at the 40 levels J = 61,62,...,100 and is given 


moet ne form 


ial 
Wid) = W(d) + & CREG(I,J)(SAT(I) - T(I) ] (7) 
I=1 
where W(J) 1s the sample mean mixing ratio at level J over 


1100 sample soundings that were used to determine 
the regression coefficients, 


CREG(I,J) is a set of predetermined regression coefficients 
calculated for each J-level from the 1100 sample 
soundings, 


SAT (I) is the saturation mixing ratio at the level I, 
where I = 1,2,...,l11 1S a Subset of the 40 
J-levels that is found to give the best results 
(see Table 2), and 


7 (I) is the mean Saturation mixing ratio at level I 
Peon lOO test, Soundings. 


Table 2. Eleven levels used for predictor levels in 
Water vapor retrieval 


i ow P(J) mb 
i 63 209.94 
2. 66. AEB ST) 
ae 70 299.01 
4, 77 412.26 
oe 82 509.93 
6. 90 699.03 
We 92 BE) Ly 
8. 94 810.25 
9. 96 870.35 
10. 98 258) 57) 
aN 100 1000 .00 


Thus the method uses an eleven level saturated mixing 
mero profile, SAT{(I), and converts it into a 40-level mi.xXing 


Mmieroeprofile, Wid), J = 61,62,...,100. The total 100-level 
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profile is then converted from the 40-level by using the cubic 
power-law profile due to Smith (1966) 


ne wG) Rep) (8) 


J<61 

The mixing ratio is then used in determining the water vapor 
transmissivity. 
C. ADJUSTMENT OF WATER VAPOR MASS IN A NON-HOMOGENEOUS 

ATMOS PHERE 

Assuming that the transmittance by water vapor, t[U(J)], 
in a spectral interval of the VTPR is a known function of the 
Quantity of the precipitable water vapor U, temperature T, 
and pressure P, an accurate approximation to t{U(J)] in the 
form cf a least squares fit polynomial can be caiculateda. 
Weinreb and Neuendorfer (1973) proposed such a best-fit model 
which in effect becomes an analytic function for calculating 
homogeneous path transmittances. In the Weinreb-Neuendorfer 
method, an inhomogeneous atmosphere is treated as a sequence 
of homogeneous layers, each having a constant temperature, 
pressure, and absorber mass. The transmittance model provides 
pera homogeneous path as a known function of T, P, and U, 
where U is the effective water vapor mass. If t(J) is the 
transmittance in the spectral interval of the instrument at 


level J, then 


t(3) = t(P(d),T(d),UI)1 . (9) 
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The procedure involved in defining t[(J)] at each level 
as depictéd in Fig. 6, uses the functional expression, Eq. 15, 
for layer 1 without modifying the precipitable water vapor 
mass, U, from its initially specified value based upon the 
homogeneous values P(1), T(1), and U(1l). This gives rise to 


the following relationship for layer 1 


ttl) = [P(1),7(),uU(1)] . (10) 


U(1) 


U (1) 


LAYER LEVEL 


A) 





m (SURFACE) 
Figure 6. Multi-level atmosphere showing P(J), T(J), and 

iG) estructunes . 
In proceeding to level 2, the transmittance through the second 
layer is affected by the new homogeneous characteristics P(2), 
T(2), and its water vapor increment AU,, and is no longer along 


meeer fe (1) ,T(1) ,U(1)] transmittance curve defined by the 
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Senditions in layer 1 (points A,-B, eng 2.07 eet As 
necessary to determine the adjusted value of water vapor mass, 
denoted V in Fig. 7, with which the radiant transmittance 
enters the new homogeneous layer 2. This is accomplished by 


~e 


backing up on curve 2 to U = V by solving 
Pelee Pe (2) 0 (2), Vv] (11) 


for V. The solution for V is accomplished by the Newton- 
Raphson iteration. 

This corrected starting water vapor mass, V, with which 
to enter layer 2, has been deduced using curve 2 in Fig. 7. 
For the purpose of computing the transmittance at level 2, the 
new total water vapor maSs is then taken as 

W = V + AU (12) 

where AU. is the known increase in precipitable water vapor in 
layer 2. This moves the transmittance value tT[P(2),T(2) ,W] 


down to point Ba on curve 2 where the value of 1(2) is then 


taken as 
t(2) = TP (2) ,T(2) ,U(2) = (13) 


In proceeding downward to level 3 through the next homogeneous 
[imwer Of Characteristics P(3), T(3), AU, the corrected amount 


ed 


of water vapor at level 2 is revised from U(2) to V which 
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T= 1 
1 
aaa 
t 
T (Pt, U) 
CURVE 2 
I 
; T( P2120) 
i 
i 
T=0 
Figure 7. Transmittance curves for atmospheric layers 1 and 2 


illustrating modification to effective U(J) ina 
non-homogeneous atmosphere. 
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becomes the level 2 starting point of water vapor mass. 


Hence the new function 


1(3) = tIP(3),7(3) ,U(3)] 


s (14) 
U(3) = V + AU 


3 
is solved for the new starting mass at the top of level 3. 
The procedure continues aS a Series of Newton-Raphson itera- 
meme tO CcOrrect the effective water vapor masses for the 
radiant transmittance upon succeeding layers. This procedure 
1s followed through each level J = 2,3,...,99 bounding the 
top of the 98 layers below J = 2, 

The procedure iS shown by Weinreb and Neuendorfer to be 


exact for a Beer's-Law monochronatic transmittance function. 


D. WATER VAPOR TRANSMITTANCE MODEL 
The water vapor tranSmittance calculations use the method 
of Weinreb and Neuendorfer (1973) applied to a modification 
of the polynomial representation of water-vapor transmission 
of Smith (1969). The basic equation for the line transmit- 
tance of water vapor implies that 
14 


eee tk = fe. xX, PGs) 
I=1 
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Here 


fl 8 eee 
= Ut i 
Xo 7 Con |} In5a4 Xo = aX) 
X4 = Jn(P/1000) X16 = XX. 
Xy iene 727.3) Xa = X 4X2 
x. = X.X epee 
5 Doe 12 4 
a 413 7 X3X_ 
nag 7 
Xa = (X5) Xa XX 


and C_. is a set of coefficients that were predetermined by 
Fitting Ec. 15 by the least-squares method using line-by-line 
calculations of the transmittances for 130 homogeneous (con- 
Stant pressure and temperature) paths. From Eq. 15 it follows 
that 


14 


t(U(J)] = expl[-exp( a Ce) (16) 
Thess 


PHoevicalucs Of X. Ehrough xX are directly calculated only 


3 14 
1f U is known. However, Xo must be calculated indirectly 

Since U, the water vapor absorber mass, iS essentially a sum 
of scaled values of AU between SuccesSSive J-levels. Newton's 


method was used to solve for the new value of X51 and subse- 


quently other Xs = X.Xo, using the first iteration of 
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() 


a se (NO U. as the first gueSs to arrive at the effective 


une? at level 2, etc. on to the subSequent J-levels. 

Water vapor transmittances in the infrared region have 
been found to contain a spectral feature not previously ex- 
Plainable by the classical line-Spectral approaches. Bignell 
(1970) found a continuum absorption from the edge of the 20 
um rotational band which overlaps the lines in the water-vapor 
and carbon-dioxide absorption regions. The continuum absorp- 
tion coefficient was found to be due to the presence of the 
water vapor dimer (H,0) >. The dimer absorption coefficient in 
the various spectral regions was found to be proportional to 
aa) (see Eqs. 17 below). Thus, it iS necesSary to add 
the dimer correction to Eq. 15 to fully account for the com- 
bined transmittances of water vapor in all channels. The 
dimer contribution at level J to the In[-lnt(J)] of Eq. 15 
due to water vapor is given by the adjusted dimer mass profile 
in terms of 


A(J-1) +A (J) 
2 


DIMER(J) = DIMER(J-1) + [ Wle(s) =P (J-1)) 


(iF) 


=" 


303.16, 


A(3) = P(3) (wi) 1° [a 


The actual transmittance due to the effective dimer mass 


in each channel is then given by 
In[-ln t(DIMER)] = B(1I)*DIMER(J) (18) 


oy) 





mere B(l), IT = 1,2,...,8 18 a set of 8 constant absorption 
coefficients, the values of which have been given for each 
channel after Bignell (1970). The combined water-vapor 


transmittance in each channel then becomes 
t{U(J)] = exp{-exp[DIMER(J)*B(I)]} t(J) (U9) 


where t(J) is the value computed using Eq. 16. 
The zenith angle correction for water vapor transmittance 
Pomincluded in DIMER(J) and U(J) of Eqs. 9 and 17 by making 


the following path-corrections for mass: 


DIMER(J) = DIMER(J) *PATH, 


U (J) Digs) * PATH; (20) 


PATH = sec @Z. 


Here sec Z is the secant of the zenith angle from the point 
representing the center of the box to the radiometer (see 


Figure 4). 


mee 6 CARBON DIOXIDE TRANSMITTANCE MODEL 

A polynomial expression for the carbon dioxide trans- 
Pmmetance 1t{(C(J)] in terms Of logarithmic functions similar 
to Eq. 15 has been shown to be feasible by Smith (1969). 
However, the absorber mass of carbon dioxide is proportional 


to pressure, and little dependence upon higher powers of the 
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in(P/1000) term has been found necessary to fit line-by-line 
computations of transmissivities with increasing pressure in 
nonhomogeneous atmospheres. 

However, the carbon dioxide polynomial expansion, similar 
to Eq. 15 for the water vapor model, has been found to be 
subject to path correction primarily for the temperature varia- 
tion along the vertical. The method of Weinreb and Neuendorfer 
(1973) was shown by Fleming (1974) to be converted into the 


following form 


ieee ice jilats)+e(n) AT (a)+y (gd) AT- (a) 1 
(21) 
+ Da he US SAL Ste iy, 


where t[{C(J)] is the transmittance due to carbon dioxide at 
level J: 


AT(J) is the temperature difference between the U. S. 
Standard Atmosphere (1962) and the average 
temperature over the layer bounded by the 
pressure levels J and J-l (J = 2,3,...,100), 


a,B,Y are a set of 3 by 100 transmittance factors that 
are tuned specifically to the filters of the 
instrument in use for the data period, 


6. is a Set of six adjustment factors applicable to 
systematic corrections for temperature varia- 
EVONSMEGOM tte Initial standard, and 


t*(J) is a correction term approximated after the 
Newton-Raphson method from 


Z 
D(J) [AT (J) ] op 
[AT (J-1) ]* 


<*(J) = A " 
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where for the initial level J =1, A, B, D are defined as 


ae, a (1) 


Baa) pa Ar (ie) (23) 


2 


p(1) = (SEL + y (a) ) (or (1) J 


and for subsequent layers, J > 1, the values of A, B, and D 


are updated iteratively by 


1) sy) = at?) (g-1) 


A 
(2) 
(1) B (J-1) AT (J) 
B (J) Peay (24) 
Ape 2 
p'1) (3) _ D (J Seis 
[AT (J=1) ] 
(2) pb’) (3) tar (gs) 120(J) (1) 2 
D (J) = +--+ A (J) y (J) [AT (J) ] 


(AT (J-1) ] 
B'2) (53) = B|) (syacay + a™ (3) B(3) arcs) 


a (2) a‘) 


(J) a(J) 


(J) 


all of which must be computed in the given order. 


Ir. OZONE TRANSMITTANCES 
MiemezoncetLansmittances t1OZ(J)}] are minor contributions 


appearing in the form of Eq. 3 for each carbon dioxide channel. 
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They are constant in each channel at a given level, varying, 
however from level to level. manus the etfect of 1T[OZ (J) } 
is expressible in the form of a 6 by 100 array, one member 
of the array for each carbon dioxide channel at each J-level. 
The channel values of tT[0Z(J)] used here are listed in Appendix 
B, as extracted from the NESS CLRAD program. 

The t{OZ(J)] values in Appendix B have not been path 
corrected and are based upon globally averaged NTP absorber 


paths of ozone in the carbon dioxide channels. 


G. TRANSMISSIVITY ADJUSTMENT FOR ZENITH ANGLE 

When the transmittance in the spectral interval of the 
VTPR has been calculated, an adjustment must be made to 
Meee. for the path differences that cccur with varicus 
zenith angles. The VIPR is designed to meaSure clear-column 
radiances from points which sSubtend a zenith angle of 23°47! 
left, 0°, or 23°47' right of the satellite track. These 
angles correspond to the center of the three boxes shown in 
Figure 3. The adjustments to the calculated transmittances 
must be made for these basic angles and any small deviations 
from them. 

The principle of the correction stems from the fact that 
if the zenity angle is other than zero degrees, the atmospheric 
absorbers will have nee effect due to the geometry involved 
(Figure 4) with the longer path, and, therefore, the trans- 


missivity in the layer will be decreased level for level. 
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The correction for both the carbon dioxide and ozone 
transmittance at level J is taken from a section of the NESS 


CLRAD program [McMillin(1974)] and is in the form of 


re = (os) >) 


A = (PATH) ea, 


(26) 
where PATH is the secant of the zenith angle, and B(J) is a 
polynomial function for carbon dioxide alone. For the ozone 
path corrected transmittances B(J) = 1 for all atmospheric 
Mevels. B(J) is a function of ln{P(J)] with terms in In[T(J) ] 
up to quadratic (and a product of these logarithmic functions), 
with Bese nO. five predetermined constants for each channel. 
These constants have been tuned for the particular VTPR 


instrument in use during the period. 


42 





ieee ey Al OF TEMPERATURE PROFILES 


PeeeeeN PUT DATA USED FOR TEMPERATURE RETRIEVAL 
1. First Guess Temperatures 

hmocotateneco the climatological soundings, T(P), 
from the U. S. Standard Atmosphere Supplements (1966) (See 
Appendix B), the retrieval program requires a first guess 
temperature profile in the box center for each VTPR clear- 
columm radiance set. The first guess temperatures are of 
two basic types depending upon the latitude of the sounding. 
For tropical latitudes, from the equator to 18 degrees North, 
radiosonde data is averaged over all latitudes and longitudes 
to provide a single climatological first guess profile. For 
the other northern hemisphere retrievals, the first guess pro- 
file comes from an application of the National Meteorological 
Center (NMC) prognostic model for a period of 12 to, at most, 
18 hours from the preceding complete analysis. The profile 
is extrapolated in pigs mad space to the VTPR clear-column 
Beemcading Location. 

The first guess temperatures, as read into the program, 
give temperatures at the 15 mandatory levels shown by an aster- 
mene (*) in Table 3, from 100 mb to 10 mb. The U. S. Standard 
Atmosphere Supplement climatology value is adjusted for latitude 
and month and then fitted to the first guess profile at the 10 
mb level. The adjustment, defined as first guess (10 mb) minus 


standard (10 mb), is then added to all values of the standard 
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Table 3. 56 level atmosphere with merged climatology and 


15 level (*) first guess field for sample sounding. 


LEVEL PRESSURE (mb) A) BEND SUE Up geyeg (G3) 
1 OvscEl S78) (en eeie 
23 OO “04.311 
34 OOO tS S100) 
35 15.00 ita), S190 
36 205. 00 = i oo 
yi ep UG, yo) ee OY 25. 
38 S00 0% = 310 . G00 
oe 40.00 Soi oo OO 
40 au. 00 = —6 30.0.0 
Al 60.00 -65.600 
42 20. 00 6S. 
43 85.00 162) 5 4a 
44 90.00 =69 20100 
45 LOG 0.0 * = EGON, 
46 SO 0.0 = =G316 400 
47 200-2002 -58 .200 
48 Zoe, 00* -48.400 
49 3007002 aoe 1.o) 7 w110)10) 
50 400.00* = 1, 
51 B00 500< ld 131048 
52 600.00 (0) 2) 3010, 
53 700°..007 4.800 
54 S000 * 13.600 
55 900500 WGaZo7 
56 OOO. 00% Z6.0:0 


44 





atmosphere, above 10 mb, and the profile is expanded to the 
56 levels shown in Table 3. 

The 56 level temperature profile is then expanded to 
the 100 J-levels, where the pressure at level J is given by 


C/2 
Peo omit (id) (0. 26087836) | 27) 


where J = 1,2,...,100 and defines the pressure from 0.01 mb 
me, 1000 mb. 
2. Other Input Data 

The remaining input data required, having given the 
VTPR first guess sounding, includes the zenith angle, sea 
surface temperature, and clear~column radiances. This data, 
along with all of the verifying radiosondes and first guess 
temperature profiles to 10 mb was kindly provided by NESS in 


the form of an archival tape. 


PeevAtTHEMATICAL DEVELOPMENT OF THE TEMPERATURE RETRIEVAL 
Given a non-scattering atmosphere in thermodynamic equili- 
brium, the clear-column spectral radiance at the top of the 
atmosphere in each channel is related to the vertical tempera- 
ture profile and the absorber mass Structure by the RTE, 
Ea. 1. 
The right hand side of Eq. 1 is composed of an earth's 
Surface radiation term, and an atmospheric contribution which 
is integrated from the surface of the earth (1000 mb) to the 


mem Of the atmosphere (0.01 mb). 
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When the 100 J-levels developed in Eq. 27 are used, Eq. l 
can be rewritten 
J=1 ay 


I, = B,[T(100}]1, (100) + fs (B,Ir(a)) _ a‘? 


: J=100 ag) a 


(28) 


The Planck radiance function used in Eq. 28 at level J and 
in channel i is 


3 
(Cl) *v; 


Cxp lasy! - l 


where ve is the wave number in channel i, 
Cl is the constant 1.9061x10 > erg em- sec + ster -, and 
G2 1s the constant 1.43868 cm °K. 

In order to reduce the effect of truncation error in the 
Mmicegration function of the RTE, and increase overall stability, 
the 100-level atmosphere was reduced to a 33-layer atmosphere 
as shown in Fig. 8. In integrating the RTE over the 100 level 
atmosphere, Simpson's rule was applied over each sSublayer K 
centered at the fractional J levels of the K-~layered atmosphere 
(Fig. 8). The main virtue of the Simpson quadrature scheme 
uSing the weighting factors of 1, 4, 1 over each layer (e.g., 
we Ol, 02.5, 04), lies in the small discretization error, 
(1.5)>/90 times the fourth derivative with respect to J of the 
Miitecrand Lunction Of Eq. 28, in substitution of the finite 


Summation for the Riemann integration. The determination of T 
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1000 MB 
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Mmegree oS. Thirty three K-layer atmosphere model used for 
the integration of the Radiative Transfer Equation. 
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mabe vels J = 02.5 etc., 1S accomplished by interpolation of 
the T(P) using the pe/7 pressure levels, defined by Eq. 27. 
Therefore, Eq. 28 can be rewritten 

3:3 


: B, (K) At, (K) : (30) 


I. = B. [T(100)] tT; (100) + ; 


K 
where 

= ae 

B,(K) = ae a l=) 4 4B. [T(K) ] ~ BE Col esl) te) 


(oak) 


is the layer-mean Planck function for layer K employing 
Simpson's rule. Moreover At ; (K) now spans three J levels 


arg. 8), 
Mp CS) =) ane (ea tRS Sy Gg (pas), (22) 


where Jy is the value of J at the center of layer K. 


Using the inversion methcd of Smith (1970) as modified by 
Martin (1974), the retrieval of mean layer temperatures is 
Getermined from clear-column radiance values as follows: 

1. Determine the first guess radiance value in 
channel i, a, sem EC. 3 Or 


2. Form the residual by comparing the clear-column 


radiance with the computed value from Eq. 30, using 


(1) _ (1) 
Pe 95 2=CLEAR ~ 72 ae 








See remenercestaual 15 not close to zero the first 
guess temperature profile ‘ae Coe ) must be improved to 
K 


ri?) (p. |) Sn aecostbvace neon 
K 


(2) E (1) (1) 
B (ree (Py V1 = Bit’ (Py 1 + Ry (34) 


This tends to move the temperature profile closer to the final 
desired convergent profile. If n steps are required to 


converge the residual close to zero, the final step is 


Geer IL) _ (n) (n) 
B. [T (P, )J = B. [T (P, yy ct R, (35) 
K K 
where finally 
RY”) -4 
~———— x 110 (36) 
i-CLEAR 


When a convergent Planck function is determined for each 
channel, i, and at — level, J, the temperature at that 
level can be. determined by solving Eq. 29 for T(J). The 
value of the Planck function in this casSe muSt be a weighted- 
mean Planck uSing the six values of B. (J), ele es, Oy ate 
each J-level. Figure 9 shows the values of transmittance in 
each of the six channels with respect to pos? Ghat is. wath 


respect to the linear J-scale. Figure 10, the first derivative 


with respect to J of the transmittances in Fig. 9, is the 
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CO, TRANSMITTANCES 
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700 
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Figure 9. Carbon dioxide transmittances channels 1 through 6, 
SE cemcMiiiam et al. (1973. 


weighting function used in determining this weighted mean 
Planck. The uSual procedure as uSed by Smith et al. (1972), 
involves the use of a constant normalized wave number 

(7 = 700 em71) imc si pewhereas in this thesis, a normalized 
wave number, decreasing with decreasing values of pressure was 
found to give a more accurate reSult, Similar to results of 
Moran (1974). By using this varying wave number (Table 4), 


advantage can be taken of the level of maximum contribution of 
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Table 4. Normalized wave numbers and corresponding pressure 
Wevelts weed im inverting the Planck function for 
temperature retrieval. 


Pressure P(J) Wave Number (v) 
mb cm 
1000 - 351 720 
350 - 96 710 
95 - 0.01 690 


CO, WEIGHTING FUNCTIONS 





2 
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LW 
x 
— 
Sp) 
1S) 
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a. 
o 2005 200 015 026 025 .030 035 .040 
d¢( TAU) /dx 
Figure 10. Weighting functions used in determining a 


weighted mean Planck function. 
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that channel to the derived equivalent blackbody temperature. 
The weighted mean Planck function then becomes 


SS) ae 
a B; (J) At, (J) 


(n) 
B [v,T(J)] = —————_____—_—___.. (37) 


- 


| 1 oy 


ao ) 
1 a 


at 
By uSing this weighted Planck, the 100 mean temperatures, 


T(J), can be determined by the inversion of Eq. 29 


T(J) = i: ne Se (38) 


mie 6as. 3/ and 38, which are accomplished after the iteration 
procedure in each channel is complete, the values of At; (J) 
are centered differences taken over AJ = 2, but centered on 
each individual J level. A forward AT. difference is used to 
initiate the process at level J = 1 and a backward difference 
is used to terminate the weighting process at J = 100. In 
mais Solution procedure, EQS. 37, 38, the algorithm used is 
essentially after Smith et al. (1972), except for the use of 


the variable v. 


C. STEPWISE TEMPERATURE RETRIEVAL 

The procedure for retrieving the layer-mean temperatures 
fmmelaver J = 1,2,...,100, from the clear-column radiances 
can be summarized as follows (program names refer to the 


Misting at the end of this thesis). 
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1. Input the necesSary data. 
a. sounding latitude. 
b. zenith angle. 
c. 15-level first guess temperatures. 
d. sea surface temperature. 
e. clear-column radiances in channels 1 to 6. 

2. Expand the 15-level first gueSs temperature 
profile to 56 levels and merge it with climatology at the 
10 mb level using subroutine UPPER. 

3. Expand the 56-level temperature profile to 100 
levels using Eq. 27 and Subroutine LEVELL. 

4. Calculate the atmospheric transmittance due to 
carbon dioxide in channels 1 through 6 at each level J uSing 
Pm@erOoucine RDTEIMP ana adjust for the zenith angle by using 
Subroutine MARTAU. 

5, Calculate the water vapor content of the atmosphere 
in the form of a mixing ratio profile, one value at each level 
Mec ing function GOFF and Eqs. 7 and 8. 

6. Usethe mixing ratios to calculate the path 
corrected atmospheric transmittance due to water vapor in 
channels 1 through 6 at each level J using subroutine TRANW. 
(Note: steps 5 and 6 were performed on the first iteration 
only, due to the small change in water vapor transmiSsivity 


from the first guess to the final temperature profile. 


a 





ieueate uate the Lotal atmospheric transmisSivity 
in each channel and at each level through the product of 
Eq. 3 which uses path corrected values of transmittance due 
to carbon dioxide, water vapor, and ozone. 

8. For the first iteration, uSing the 100-level 
first guess temperatures, the Planck functions must be 
calculated in each of the 6 channels. 

a. Calculate the 100 J-level Planck functions 
meng Bq. 29. 

b. Calculate the 33 K-level Planck functions by 
fitting a fourth order Lagrangian polynomial to the 100 J-level 
Planck-profile [after Conte and deBoor (1965)]. 

9. Calculate the At(J) for each channel in the 100 
J-level and At(K) for the 32 K-level atmospheres. The Az(K} 
Will be used in the discretized form (Eq. 30) of the RTE, in 
Order to reduce the truncation error associated with the 
finite Summation, Substituted for the Riemann integration of 
fmemimetegrand function in Eq. 28. This is important while 
dealing with the t(J) values which occur before the temperature 
profile and its associated transmisSivities have converged. 
The converged At(J) will be used with Eq. 37 in producing the 
new 100-level temperature profile. 

10. Compute the layer-mean Planck values by Eq. 3l. 

11. Compute the radiance in each channel uSing Eq. 30 
and determine the residual difference between the Bienes colin 


radiance and the computed radiance, Eq. 33. 
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2. Lf the residual value is not close to zero, 
another iteration must be performed to improve the accuracy 
of the retrieved temperature profile. 
a. Calculate the new 100-level and 33-level 
Planck functions by adding the residual to each previous 
value, Eq. 34. 
b. Calculate the new 100-level temperature 
Peerate using Eqs. 37 and 38, and go back to step 4. Continue 
recycling until convergence is achieved (see Eq. 36). 
13. If the residual is small enough (Eq. 36 is true), 
go directly to the calculation of the new 100-level temperature 


eeorile, Eqs. 37 and 38. 


pee ROGRAM OUTPUT 

The program is designed to produce a mean layer temperature 
profile at 100 levels from 1000 mb to 0.01 mb. For the purpose 
of testing the accuracy of these retrieved temperatures, 
additional inputs were needed. 

First a set of radiosonde soundings, as close as possible 
to each retrieval site in latitude, longitude, and real time 
was read in. While these soundings are not considered to be 
the exact state of the atmosphere at the retrieval site, they 
are as close as possible and therefore were used to determine 


the error in the retrieval program at each level J by using 


ERROR(J) = T(J) - AT(J) (29) 
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where T(J) is the retrieved temperature, and 
AT(J) is the radiosonde sounding temperature. 

Secondly, a set of retrieved temperatures Supplied from 
NESS for each case investigated was read in and compared with 
the radiosonde temperatures. 

The program then output the following for each of the 
233 retrievals attempted (see sample output following the 
Program listing): 

meeotation number (1,2,...,233) 

2. Sounding zenith angle in degrees. 

See Station latitude. 

4, The adjustment at 10 mb necesSary to fit the first 
guess temperature profile to the climatological 
temperatures. 

5. The number of iterations necessary for convergence. 


6. Table of results, listing for each of the 15 levels: 


a. Pressure 


b. Retrieved temperature 
c. Radiosonde sounding temperature 
d. Error (retrieved -- sounding) 


e. NESS retrieved temperature 
f. Error (NESS retrieved - sounding) 


In order to easily compare the accuracy of the retrieval 
with the radiosonde data and also to compare the error of 
this method with that of NESS, histograms were drawn and 
Statistical data were obtained. The errors were compared 
level for level on an overall basis (globally averaged) and 
also on a latitude-band basis for the latitude ranges shown 


Mae Table 5, 


56 





Table 5. Latitude bands used for error analysis of the 
VTPR retrieved temperatures. 


LATITUDE BAND BAND WIDTH 
iwi 0 < LAT < 20N 
2 20N < LAT < 40N 
3 40N < LAT < 60N 
4 60N < LAT 


7, 





Veeeekecoublo OF THE VITPR TEMPERATURE RETRIEVAL 


In general, the errors (retrieved temperature - radiosonde 
temperature) were less than 5 °K when averaged over the 
latitude-bands shown in Table 5. The errors also compared 
favorably with the results of the NESS VTPR retrieval in all 
15 atmospheric levels and the four latitude bands. In 
calculating the errors, the assumption that the radiosondes 
matched the VTPR retrieval sites exactly in time and space 
waS made. The difference, though small (averaging +0.6 
degrees in latitude, +0.7 degrees in longitude, and +1.87 
hours in time), coupled with the small data base of 233 
retrievals, could have had some effect on the results and 


Smould be of some consideration. 


Ae GLOBAL RESULTS 

The global mean results obtained here are summarized Tea 
Table 6 along with the averaged results from the NESS 
retrievals. The frequency distributions of the sileseil error 
analysis at each of 15 key atmospheric levels are shown as a 
function of the error classes encountered (See Appendix A-1 
to A-15). These frequency distributions, in the form of 
histograms, are included as a visual aid in analyzing the 
deviation of the error from the mean. They show that even 
though the range of errors can extend over as much as 20 °K, 
the probability of a relatively large error is much less than 


the probability of the small errors in the vicinity of the mean. 
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Table 6. Retrieved temperature error analysis on a global baSis., 


RETRIEVAL NESS 
PRESSURE (mb) MEAN RMS ERROR MEAN RMS ERROR NO. OF 
ERROR (°K) ERROR (°K) SAMPLES 
10 -0.55 4.00 -0.75 4 Tl 59 
20 0.57 3.25 -0.11 Pe 2 ey 
30 -0.06 2.20 -0.17 Be 12 42 
50 Oo7 3.20 -0.17 ol 142 
70 0. WE 4.39 0.91 5 ue 179 
100 -0.45 3.89 -2.30 5.12 179 
150 -0.30 2S -1.51 CW ke 208 
200 -0.66 3.33 -1.04 3.29 215 
250 215 ee 0.41 3.33 206 
300 -0.58 Deo 0.24 27/0 223 
400 -0.18 219 -0.09 2.07 226 
500 -0.67 3.09 -0.32 2.13 AG) 
700 -0.99 3.29 -0.86 2.81 226 
850 -0.08 BG 2E -0.31 go> 227 
1000 -0.24 eed 0.26 3.89 196 
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Bee LATITUDE BAND RESULTS 

Because of space limitations, frequency distributions 
were not included for the latitude band results. However, 
Table 7 furnishes a Summary of the latitude band analysis 
of errors including their root mean square errors along with 


the corresponding NESS results. 
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VI. CONCLUSIONS 


The first objective of this thesis was to incorporate the 
State-of-the-art methods [as used by Fritz et al. (1972)] in 
calculating atmospheric transmittances, into the retrieval 
Scheme. The method involved here for calculating the effects 
of water vapor, ozone, and carbon dioxide transmittance proved 
to be satisfactory. This, coupled with the inclusion of a 
33 level integration scheme and the use of a variable wave 
number, aS opposed to the constant v = 700 ai, combined to 
Give reSultsS comparable to those of NESS. 

The temperature profile convergence occurred after an 
average of 5 iterations. While this is an improvement over 
the Moran (1974) model, as would be expected with the more 
sophisticated procedures incorporated, the number of iterations 
required would possibly be reduced further through the use of 
a more accurate first guess temperature profile. Gelman et 
al. (1972) used a regression technique for estimating the 
first guess temperatures above 10 mb, aS opposed to merging 
with the standard atmosphere climatology as used here. Since 
this regression technique of Gelman produces a more accurate 
first guess temperature profile, it would reduce the number 
of iterations required for convergence even further. However 
this program was not available locally for inclusion in this 


emesis .« 
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The second objective, that of evaluation through compari- 
son with matching radiosondes, showed that the retrieval 
system produced errors that were comparable to NESS retrievals, 
Or better, especially in the area of the tropopause. The 
small data base of 233 soundings is not considered large 
enough to evaluate the accuracy of the retrieval system beyond 
that, and more experimentation with a larger sample, and 
inclusion of the Gelman regression technique is required. 

It is concluded that this thesis has demonstrated a 
consistent and accurate means of temperature profile retrieval 


from VTPR clear-column radiance data. 
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APPENDIX A-l 


Histogram’of the globally averaged errors showin Lika 

| g probabilit 
emstribution [P(E) ] of the occurrence of the Curors,. the eae, 
error frequencies in the error class E (marked along the top) 
the Mean (M), and the Quartiles (-), for sample size 196. | 
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APPENDI*X A-2 


Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the Errors, tne actual 
error frequencies in the error class E (marked along Ene top) - 
the Mean (M), and the Quartiles (-), for sample size 227. 
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APPENDIX A--3 


Histogram of the globally averaged errors showing probability 
emseribution [P(E) ] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along Ene top) ,; 


the Mean (M), 
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APPENDIX A-4 


Histogram of the globally averaged errors showing probability 
Petri pbution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles (-), for sample size 227. 
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APPENDIX A-5 


Histogram of the globally averaged errors showing PEObaen li) ty 
distribution [P(E)] of the occurrence of the GE EOrmsS ache aceua | 
error frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles (+), for sample size 226. 
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APPENDIX A-6 


Histogram of the giobally averaged errors showing probability 
Mme oution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles (-), for sample size 223. 
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APPENDIX A=-7 


Histogram of the globally averaged errors showing probability 
distribution [P(E)]} of the occurrence of the errors, the actual 
meeeor <trequencies in the error class E (marked along the top), 
BeeeMean (M), and the Quartiles (-), for sample size 206. 
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APPENDIX A-8 


Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the CErers, thewactual 
error frequencies in the error class E (marked along the top) 
the Mean (M), and the Quartiles (-*), for sample size 215. 
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APPENDIX A-9 


Histogram of the globally averaged errors showing probability 
Seeetribution [P(E)] of the occurrence of the errors, the actual 
error frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles (-), for sample size 208. 
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APPENDIX A-10 


Histogram of the globally averaged errors showing probability 
distribution [P(E)] of the occurrence of the errors, the actual 
SyeeOr frequencies in the error class E (marked along the top), 
Ene Mean (M), and the Quartiles (+), for sample size 179. 
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APPENDIX A=-11] 


Histogram of the globally averaged errors showing probability 
@istribution [P(E)] of the occurrence of the errors, the actual 
Serer (frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles (-), for sample size 179. 
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APPENDIX A-12 


Histogram of the globally averaged errors showin bala 

. G prebacilit 
erstribution [P(E) ] of the Sceurrence OF the errors, the Sai 
error frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles ae tOr Sample size 42. 
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APPENDIX A-13 


Histogram of the globally averaged errors showin L1i 

2G eigleloevonbilitie 
distribution (P(E) ] of the Occurrence of the CrrOne .aene Been 
error frequencies in the error class F (marked along the top) 
the Mean (M), and the Quartiles (-), £or sample size 142. 
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APPENDIX A-14 


Histogram of the globally averaged errors Showing probability 
distribution [P(E)] of the occurrence of the errors; the actual. 
error frequencies in the error class E (marked along the top), 
the Mean (M), and the Quartiles (-), for Sample size 112. 
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APPENDIX A-15 


Histogram of the globally averaged errors Showing probability 
eeestribution (P(E) ] Of the occurrence of the CrrOrs tie actual 
error frequencies in the error class F (marked along the top) 
the Mean (M), and the Quartiles (-), for sample size 59. 
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